Majorana fermions in s-wave noncentrosymmetric superconductor with Rashba and 

Dresselhaus (110) spin-orbit couplings 



Jiabin You 1 ^ C. H. Oh 1 ^ and Vlatko Vedral 1 ' 2 ^] 

1 Centre for Quantum Technologies and Department of Physics, 
National University of Singapore, 117543, Singapore 
2 Department of Physics, University of Oxford, Clarendon Laboratory, Oxford, 0X1 3PU, United Kingdom 

(Dated: September 6, 2012) 

The asymmetric spin-orbit (SO) interactions play a crucial role in realizing topological phases in 
noncentrosymmetric superconductor (NCS). We investigate the edge states and the vortex core states 
in s-wave NCS with Rashba and Dresselhaus (110) SO couplings by both numerical and analytical 
methods. In particular, we demonstrate that there exists a novel semimetal phase characterized by 
the flat Andreev bound states in the phase diagram of the s-wave Dresselhaus NCS which supports 
the emergence of Majorana fermion (MF). The fiat dispersion implies a peak in the density of states 
which has a clear experimental signature in the tunneling conductance measurements and the MFs 
proposed here should be experimentally detectable. 
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Topological phase of condensed matter systems is the 
quantum many-body states with nontrivial momentum 
or real space topology in the Hilbert spaces. A list of 
the candidate systems includes superfluid 3 He pQ, quan- 
tum Hall systems [H [3] , topological insulators (TI) [U- 
H], topological superconductors (TSC) [3 [TU] and so on. 
Theoretical predictions and experimental observations of 
TI and TSC have spawned considerable interests recently 
as a promising tool for spintronics and topological quan- 
tum computation (TQC). The TIs are insulators on the 
inside and conductors on the surface, occurring in the 
presence of spin-orbit (SO) interactions and time-reversal 
symmetry (TRS). The SO coupling locks the momen- 
tum and spin of the surface electrons and the TRS com- 
pletely cancels the surface charge current and leads to the 
pure spin current, as such the TIs open a new gateway 
for spintronic devices [TlTU3| . Another peculiar prop- 
erty of the topological phase is that it hosts Majorana 
fermion (MF) [9j [TO] HH [15] which is its own antiparti- 
cle. Discovery of these MFs would be an achievement 
in itself, and could also lead to a potential breakthrough 
for technological advance in the fault-tolerant TQC [T5] 
which may overcome the decoherence, the main obsta- 
cle to realize a working quantum computer. There are 
several proposals for supporting MFs in condensed mat- 
ter systems, for example, chiral p-wave superconductor 
[§] , superconductor-topological insulator-superconductor 
(S-TI-S) junction [TS], TSC QH \U\ EB] and so on. The 
TSCs have been first experimentally realized in physical 
systems of the ID superconducting indium antimonide 
nanowire [B3 and the Andreev bound states (ABSs) of 
Cu :E Bi 2 Se 3 gOlllI]; the signatures of MFs in S-TI-S struc- 
ture has also been reported in recent experiment |22j . 
The S-TI-S junction and TSC could pave the way for 
the implementation of the TQC and pending to be fully 
explored theoretically and experimentally. 

In this paper, we study the MFs in the edge states and 



the vortex core states of a kind of TSC named s-wave 
SO coupled noncentrosymmetric superconductor (NCS). 
Here we focus on the Rashba and Dresselhaus (110) types 
of SO couplings which are induced by the structure inver- 
sion asymmetry and bulk inversion asymmetry, respec- 
tively. It is found that the asymmetric SO interaction 
plays a crucial role in realizing topological phases in the 
NCS. Although the Rashba SO coupled NCS has been 
previously investigated [TU], the Dresselhaus (110) SO 
coupled NCS remains unexplored. To obtain a compre- 
hensive understanding of the impact of the asymmetric 
SO couplings on the topological phases, we mainly ex- 
plore the Dresselhaus (110) SO coupled NCS and then 
consider the interplay between them. Interestingly, we 
find that there is a novel semimetal phase in the Dres- 
selhaus NCS, where the energy gap closes at some points 
in the first Brillouin zone (BZ) and different kinds of 
flat ABSs emerge. We demonstrate that these flat ABSs 
support the emergence of the MFs analytically and nu- 
merically. It is known that the Chern number is not a 
well-defined topological invariant in the gapless region, 
however, we find that the different kinds of flat ABSs 
can still be distinguished by the number of gap-closing 
points in the first BZ. These dispersionless ABSs featured 
by vanishing group velocity can appear at certain energy, 
depending on the region of phase diagram in which we 
choose the parameters. Previously, several authors pro- 
posed the flat ABSs in NCS Li 2 Pd :E Pt3_ :E B with high 
order SO couplings [23j [24], d xy -w&ve superconductor, 
p x -wave superconductor and d xy + p-wave superconduc- 
tor [25] ■ Instead we study the flat ABSs in the s-wave 
Dresselhaus SO coupled NCS which may be experimen- 
tally more feasible. The flat dispersion implies a peak in 
the density of states (DOS) which is clearly visible and 
has an experimental signature in the tunneling conduc- 
tance measurements [22 [231 12H] ■ The zero-bias conduc- 
tance peak has been observed in recent experiments on 
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the indium antimonide nanowire and Cua;Bi2Se3 [T§1 |2"U] 
and argued to be due to the ABS. Thus the flat ABS 
and the zero-bias conductance peak in the DOS predicted 
here are experimentally accessible so that the MFs in the 
Dresselhaus NCS should be detectable. 

We begin with modeling the Hamiltonian in a square 
lattice for the two dimensional s-wave NCS with Rashba 
and Dresselhaus (110) SO interactions in a magnetic field, 
which is given by H = H kin + H z + H u + H^ 10 + H s : 

is 0—x,y is 



+ i(4-u c * - 4+u c ^) + H - c -]' 



no 
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H s =J2[A s 4 t 4 l+ R.c.], 



(1) 



where cj s (cj s ) denotes the creation (annihilation) opera- 
tor of the electron with spin s = (t, I) at site i — (i x ,i y ). 
i/kin is the hopping term with hopping amplitude t and 
chemical potential — fi. Hz is the Zeeman field induced by 
the magnetic field with components V = (V x , V y ,V z ) — 
^{B X ,B V ,B Z ). H R and H^ 10 are the Rashba and Dres- 
selhaus (110) SO couplings and H s is the s-wave su- 
perconducting term with gap function A s . We assume 
t > throughout this paper. In the momentum space, 
the Hamiltonian is H = | Ek '/'k^C 4 )^ with i/i^ = 



( c kT> c H' c -kT, c -ki), 



when 



k G 1BZ and the Bogoliubov-de Gennes (BdG) Hamilto- 



W(k) = 



& + (£k + V) • a iA s a y 

-iA s <T y -£ k + (£ k - V) • a* 



(2) 



cos k y ) — fx, £k 



{(T x ,a y ,a z ) are 



where £k = —2t(cosk x 
(a sin k y , — a sin k x , (3 sin k x ) and a 
Pauli matrices. 

The nontrivial topological order in the SO coupled 
NCS is characterized by the existence of gapless chiral 
edge state and Majorana fermion. The former is deter- 
mined by the nonzero Chern number and the later can 
be found at the edge or in the vortex of the system. We 
shall demonstrate these features in the Hamiltonian Eq. 
with different parameters. Particularly, we are inter- 
ested in the following four cases: (i) Rashba NCS in a 
perpendicular magnetic field, /? = and V x = V y = 0; 
(ii) a = P and V x = V y = 0; (iii) a = and V z = 0; (iv) 



Dresselhaus NCS in an in-plane magnetic field, a = 
and V z = 0. 

To find the topological phases in the Hamiltonian Eq. 
pi), we can use the Chern number to characterize the 
nontrivial momentum space topology of the Hamiltonian. 
The Chern number is defined as C — ^ J T2 dk x dk y F{)i), 
where J"(k) = e^'cJ^A^k) is the strength of the gauge 
field Ai(k) = iE occ .( ? /'n( k )|5/c i V'™(k)), where ip n (k) 
is the eigenstates of the Hamiltonian Eq. ([2]). The 
integral is carried out in the first Brillouin zone and 
the summation is carried out for the occupied states 
of the Hamiltonian Eq. A particular feature of 

the topological materials is their robustness against 
perturbations. As long as the topological quantum tran- 
sition does not happen, the Chern number will remain 
unchanged. Since the topological quantum transition 
happens when the energy gap closes, we can depict the 
phase boundary by studying the gap-closing condition of 
the Hamiltonian Eq. ([2| first. By diagonalizing the BdG 
Hamiltonian Eq. ([2| in the periodic boundary conditions 
in the x and y directions, we can obtain simple analytical 
solutions for the case (i) and (iv) mentioned in pervious 
paragraph. The energy spectra of them are E(k) = 

where V — V z for case (i) and V = \jv x + Vy f° r case 

(iv). Therefore, we can find that the energy gap will 
close at & + Cl + V 2 + A 2 S = 2^^ + + A*), 

which leads to the following gap-closing condition after 
some straightforward calculations, 



£ + A 2 s = V 2 ,U = 0. 



(3) 



For the case (i), the gap closes at (k x ,k y ) = 
(0, 0), (0, tt), (tt, 0), (tt, 7t). By substituting these values 
into Eq. we obtain three boundary lines, (u ± At) 2 + 
A 2 = V 2 and p 2 + A 2 = V 2 shown in Fig. fy). Inter- 
estingly, for the case (iv), by the same reason, the gap- 
closing happens in the line k x — or k x — n. Substitute 
these values into the Eq. we have two equations, 

{fi ±2t + 2t cos k y ) 2 + A 2 = V 2 , where V = ^jv 2 + V 2 . 

Associating with | cos /c y | ^ 1, we can find that the gap 
closes in the regions from A to G as shown in the Fig. 
(|T]b) . We also consider the interplay of the Rashba and 
Dresselhaus (110) SO couplings in the cases (ii) and (iii). 
There are not simple analytical gap-closing conditions for 
them. However, as long as the z component of the mag- 
netic field vanishes, the phase diagram is topologically 
identical. Therefore, the phase diagram of the case (ii) 
is topologically the same as Fig. ([T^i) and the case (iii) 
is the same as Fig. ([TJd). In the Fig. ^f), the Chern 
number is attached for different regions. We can see that 
in the case (i) and (ii), the corresponding regions have 
the same Chern number as calculated by M. Sato |10) . 
Interestingly, different from the case (i) and (ii), in which 
the gap closes in some boundary lines, in the case (iii) 



(a) (b) 

FIG. 1. (color online). The phase diagrams of s-wave (a) 
Rashba and (b) Dresselhaus NCS. 

and (iv), the gap is closed in the areas from A to G as 
shown in the Fig. ([lja), which means that the system is 
in the semimetal phase in these regions. Inside the gap- 
less regions, it is well known that the Chern number is 
not well-defined. In the following, we shall use the gap- 
closing points in the two lines k x — and k x — n to clas- 
sify the semimetal phase in Fig. ([TJa). To demonstrate 
the topological features of the s-wave Rashba and Dres- 
selhaus SO coupled NCS, we study the ABSs and MFs at 
the edge and in the vortex core of them. In the first two 
cases, the phase diagrams are topologically equivalent to 
the case (i), the Rashba NCS, which is already studied 
by M. Sato [TU1H7]. ^ n * ne l as * two cases, the phase di- 
agrams are topologically equivalent to the case (iv), the 
Dresselhaus NCS. We shall focus on the Dresselhaus NCS 
in the following paper and see what happens inside the 
semimetal phase. 

We now turn to study the edge states of the Dres- 
selhaus NCS. By setting the boundary condition of x 
direction to be open and y to be periodic, we diagonal- 
ize the Hamiltonian Eq. ([2| in the cylindrical symmetry 
and get the edge spectra of the Hamiltonian. The results 
are depicted in the Fig. (2). Interestingly, although the 
gap closes in the semimetal phase from A to G as shown 
in the Fig. (JIJd), there exist dispersionless ABSs at the 
edge of the system as shown in Fig. Later we shall 
demonstrate that these flat ABSs are MFs. The number 
and position of the edge states depend on the number 
of the gap-closing points of the Dresselhaus NCS dis- 
cussed above. From the gap-closing condition, we have 

cosfcj, = ± ^ V ~2t A -± ~ ~ ^ 1- We can find that different 
number of solutions to k y in the different regions of the 
phase diagram Fig. ^p)- If we define the number of 
gap-closing points situated in the line k x = as v\ and 
in the line k x — n as z/ 2 , then we can use {v\ 1 v?) to distin- 
guish different regions in the phase diagram. After some 
straightforward calculations, we get that (2, 0) for the re- 
gion A and C, (0, 2) for the region B and D, (2, 2) for the 
region E, (4,0) for the region F, (0,4) for the region G. 
Consequently, the edge spectrum of region A is topolog- 
ically equivalent to the region B, C and D, and similar 
relation for region F to region G. The edge spectra for 




FIG. 2. (color online). The edge spectra of s-wave NCS with 
Rashba and Dresselhaus SO interactions. The open edges are 
at i x — and i x — 50, k y denotes the momentum in the y- 
direction and k y £ (— 7r,7r]. For (a)-(c), which correspond to 
the three regions A, E and F in the Fig. ([ljj) , the parameters 
axe t = 1, V z = 0, a = 0, /3 = 1, A s = 1. The chemical 
potential and in-plane magnetic field are (a) fi = —4, V 2 = 5, 
(b) ^ = 0, V 2 = 5, (c) fi = -2, V 2 = 3. For (d), t = 1, V z = 0, 
a = p = 1, A s = 1, n = -4, V 2 = 5. 

the region A to D are shown in the Fig. ([2^), for the re- 
gion E is shown in Fig. ^p) and for the region F and G 
are shown in Fig. (|2j;). We observe that there exist dif- 
ferent kinds of the zero energy flat ABSs in the different 
regions of the semimetal phase. Especially, apart from 
the zero energy flat ABSs, we also find that there are flat 
ABSs at some nonzero energy in the region F and G. All 
of these flat ABSs have clear experimental signature in 
the tunneling conductance measurements, therefore, the 
MFs emerging in the Dresselhaus NCS should be experi- 
mentally observable. We also show the edge spectrum of 
the case (iii) in the region A and compare with Fig. ([2^,) . 
Interestingly, from Fig. ([2ji), we find that the group ve- 
locity (dE/dk) of the MFs is tunable by interplaying the 
Rashba and Dresselhaus SO interactions. 

The existence of the edge states implies the nontriv- 
ial momentum space topology in the Dresselhaus NCS 
so that the MFs will emerge at the edge of the sys- 
tem [TJ3- In the following, we explicitly calculate the 
Majorana zero modes at the edge of the Dresselhaus 
NCS in the cylindrical symmetry. Let x direction to be 
open boundary and y to be periodic, then k x — > —id x , 
we solve the Schrodinger equation of Hamiltonian Eq. 
Q in the real space, H{k x — > ~id x ,k y )^! = 0, where 
^ = (u^, it)., v-f , vi) T . Due to the particle- hole symmetry 
in the Dresselhaus NCS, we have = i?| and = u| 
at zero energy. Thus, we only need to consider the upper 
block of Hamiltonian Eq. ^ . For simplicity, we consider 
the low energy theory at k x = 0, up to the first order, we 
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have 



(s(k y ) - i/3d x )u t + (14 - iV y ) Ui + A s ul = 0, 
(e{k y ) + ifid x )u± + (V x + iV y )u t - A s u| = 0, 

where e(k y ) = — 2t(l + cosk y ) — /i. Observe that u-j- 
±ra? in Eq. (HI), we have 



(4) 



(-ij3d x + e{k y ) t iA s )u t ± (V y + iV x )u* = 0. (5) 

Solving this equation directly, we obtain that when u-j- = 
iuf , the solution is = c\u].(x) + C2M^(x), where C\ 

and C2 are real numbers and 

u\{x) = A ie XlX + A 2 e X2X , 



w^(a;) = iBie Aia; +iB 2 e A2;l: , 



(6) 



Vx-i(V B +e) - 



2_c-2 V' 



where Ai j2 = ~ A °TV^-^ and A 1)2 = 

#1,2 = |(1 ± ^l^' ) i when u-f = — the solution 
is u-f(a;) = C3U^(x) + C4U*(x), where C3 and C4 are real 
numbers and 



v%(x) = iA x e x * x + iA 2 e x * x , 
ut{x) = B ie X3X + B 2 e x * x , 



(7) 



where A3, 4 = As=FV g /2 ~ ~- We consider the Dresselhaus 
NCS in the positive x plane with the edge located at 
x = 0. Then the critical point for existing of the normal- 
izable wavefunctions under this boundary condition is de- 
termined by V 2 — s(ky) 2 = A 2 , which is consistent with 
the gap-closing condition (fi + 2t + 2t cos k y ) 2 + A 2 = V 2 
at k x = 0. By the same reason, the condition for normal- 
izable wavefunctions is consistent with the gap-closing 
condition (/.i — 2t + 2t cos k y ) 2 + A 2 — V 2 at k x = ir if we 
consider the low energy theory at k x = -k. Let's assume 
A s > for simplicity, then the Majorana bound state is 
(u^, iv,p Up ~iu^) T , where is the solution of Eq. won. 

To further study the MFs in the Dresselhaus NUS, 
we consider the zero energy vortex core states by solv- 
ing the BdG equation for the superconducting order 
parameter of a single vortex A(r, 9) = Aexp(i#) [25] . 
To do this, the s-wave superconducting term in Hamil- 
tonian Eq. ([I]) is modified to be position-dependent, 



H.c). We numerically solve the 



Schrodinger equation = for the Hamiltonian in 
Eq. where ^ = (u^, u±, v-f, v^) T . At zero energy we 
have = i>| and = v*^ as the particle-hole symmetry 
in the Dresselhaus NCS, then the Bogoliubov quasiparti- 
cle operator, ^ (E) = I^(ujt4f + U 44; + ^ C 4 + W 4 C 4) 
becomes Majorana operator 7^(0) = 7(0). Below we only 
consider the zero energy vortex core states for discussing 
the MFs in the vortex core. Let's set the x and y di- 
rections to be open boundary, then we solve the BdG 
equations numerically and calculate the density profile of 
quasiparticle for the zero energy vortex core states. The 
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FIG. 3. (color online). The probability distribution and den- 
sity of states of quasiparticle for the s-wave Dresselhaus NCS 
plotted on the 81 x 81 square lattice. The parameters are 
t = 1, V z = 0, a = 0, P = 1, A s = 1. The chemical potential 
and in-plane magnetic field are (a) /1 = —4, V 2 — 5 and (b) 
H = 0, V 2 — 5; (c) and (d) are the density of states for (a) 
and (b), respectively. 



zero energy vortex core states of Rashba NCS has been 
studied by M. Sato [TU1 [SB] - For fully understanding the 
impact of the asymmetric SO couplings on the NCS, we 
explore the Dresselhaus NCS in the present work. Pre- 
viously, we have shown in Fig. ^ that there is a novel 
semimetal phase in the Dresselhaus NCS where the zero 
energy flat edge states host MFs. Here we shall ascer- 
tain whether there exist zero energy vortex core states 
hosting MFs in this semimetal phase. The density pro- 
files of quasiparticle of the zero energy vortex core states 
are shown in Fig. (|3^i) and Fig. ^ijp), which correspond 
to the region A and E in the phase diagram Fig. ^p), 
respectively. The numerical results of the energy for the 
zero energy vortex core states are E = 2.27 x 1CP 4 for 
Fig. ^) and E = 7.11 X 1(T 6 for Fig. respec- 
tively, which are within the simulation tolerance. For 
the choice of parameters in our simulations, the order 
of mag nitude of A 2 /E F is 1CT 1 . Thus, these numeri- 
cal results have much smaller energy than the Caroli-de 
Gennes-Matricon (CdGM) mode [29]. The correspond- 
ing DOS are also shown in Fig. |3j:) and Fig. |3ji). We 
find that there is a peak at zero energy which further 
confirms the existence of MFs in the vortex core in the 
semimetal phase. This peak is measurable in the tunnel- 
ing conductance experiments and we argue that the MFs 
predicted in the vortex core may also be detectable. 
In summary, we have investigated the topological 
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phases in s-wave Rashba and Dresselhaus (110) SO cou- 
pled NCS. We study phase diagrams, edge states and vor- 
tex core states in different combination of SO couplings. 
Particularly, we find that there is a novel semimetal phase 
appearing in the phase diagram of the Dresselhaus NCS. 
The Chern number is not a well-defined topological in- 
variant for the gapless semimetal phase which however 
can still be divided into different regions by the num- 
ber of gap-closing points situated in the line of k x — 
and k x = w. We observe that there exist flat ABSs 
in the semimetal phase. The flat dispersion leads to a 
peak at zero energy in the DOS which has a clear ex- 
perimental signature in the tunneling conductance mea- 
surements. We argue that the zero energy flat ABSs at 
the edge and in the vortex core of the Dresselhaus NCS 
are MFs by analytical solutions and numerical simula- 
tions. Therefore, compare with the proposals for MFs 
in NCS Li 2 Pd I Pt3_, ! ;B with high order SO couplings, 
dxy-w&ve superconductor, p x -wave superconductor and 
d xy + p-wave superconductor, the MFs emerged in the 
Dresselhaus NCS may be more experimentally accessi- 
ble. 
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